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Abstract 

It is shown that the ratio of the production rates of the pairs B + B~ and B®B° 
should experience a substantial and rapid variation with energy within the width of 
the T(4S) resonance, crossing the value of one near the center of the resonance. This 
behavior is due to an interference of the rapidly changing with energy Breit-Wigner 
phase with the phase introduced in the wave function of charged mesons by their 
Coulomb interaction. 



The ratio of the production rates of charged and neutral B mesons in e + e annihilation 
at the T(4S) resonance, 



a(e+e~ -> B°B°) 



is an important parameter in detailed studies of the properties of B mesons. Recent dedicated 
measurements^ 12 0] of R c ^ n at the maximum of the resonance report values ranging from 
1.04 ± 0.07 ± 0.04 P to 1.10 ± 0.06 ± 0.05 2J, which leave enough room for further studies 
of the quantity of interest 8R c / n . 

Theoretically the difference 5R c / n of the discussed ratio from one arises as dominantly an 
effect of the Coulomb interaction, clearly different for charged and neutral B mesons, since 
the mass difference m^o — rriB+ = 0.33 ± O^SMeV^jl] is quite small, and its effect can be 
accounted separately. In the most simple approach j^j, where the B mesons are treated as 
point particles, and the existence of the resonant interaction is ignored, the estimate of 8R c ^ n 
can be expressed in terms of the cm. velocity v of produced B mesons, using the textbook 
Coulomb wave functions: 

At the excitation energy of the T(AS) resonance, E = My(4S) — 2mB ~ 20^^61^, one 
has v fa 0.06, and the simple estimate (J2J) would yield 5R Cj/n ~ 0.19. It was subsequently 
argued|H] that the estimate of bR c l n can in fact be substantially reduced from eq.(j2J) if one 
accounts for a finite size of the B mesons (through their electromagnetic form factor) and 
also for the finite size of the T(4S). The latter effect was further discussed in a specific 
model of heavy quarkoniump]. Recently the problem of calculation of the ratio R c l n was 
revisited jH] in the context of a chiral-type model for strong interaction of B mesons at short 
distances, including the B*Bir vertex and the coupled channels with pairs of pseudoscalar 
and/or vector mesons, although still considering all the mesons as point-like with respect to 
the Coulomb interaction. 

It should be noted that in all previous theoretical studies of the ratio R c ^ n the presence 
of the resonance in the wave function of the B meson pair was essentially ignored. In other 
words, the coupling of the resonance to the B mesons was either treated perturbatively 
(although with a form factor [HI Ej), or the considered model of the strong interaction did 
not contain a resonance at alljH]. For this reason the results predicted a smooth behavior 
of 5R c ' n with energy in the region of the T(45) resonance. It is however well known (see 
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e.g. in the two last chapters of the textbook jU]) that the presence of a resonance produces 
a large effect on the wave function of the scattering states, which is rapidly changing across 
the resonance with energy. The scale for the variation is set by the resonance width T. In 
particular, the relative phase 28 between the outgoing and incoming spherical waves (twice 
the scattering phase 8) changes by 2tt at the scale T when the excitation energy E passes 
the central value E . It is the purpose of the present paper to properly take into account 
the resonant behavior of the wave function of the scattering states along the lines of the 
standard non relativistic scattering theory [0]. It will be shown that an interplay between the 
rapidly changing relative phase of incoming/ outgoing wave and of the effects of the Coulomb 
interaction gives rise to a rather non-trivial behavior of 5R c ^ n with the energy changing across 
the resonance. Namely bR c l n has to change sign at energy within a fraction of the width Y 
from the 'nominal' resonance center energy E . It will also be argued that the effect of a 
rapid variation of 8R c ^ n should be model independent, while the details, such as the overall 
magnitude of 5R c ^ n and the precise position of its zero(s), do depend on yet unknown details 
of the strong and electromagnetic interactions of the B mesons at short distances and of 
the structure of the T(AS) resonance. Thus a detailed experimental study of the behavior 
of 5R°/ n in the resonance region could in principle provide a certain insight into those finer 
properties of the hadron dynamics. 

The standard physical picture for considering the scattering in the resonance region 
(c.f. Ref. [9 ) is that the strong interaction, responsible for the existence of the resonance 
has a short range a, and the essential effects of the interaction at distances r < a can be 
parameterized in terms of phenomenologically measurable parameters of the resonance, the 
most important being its energy Eq above the threshold and the width T. At distances larger 
than a the motion is described by a known potential V(r): either V(r) = 0, or a Coulomb 
potential (with a possible modification due to form factor at short distances), where the 
wave function of the scattering state can be found explicitly from the Schrodinger equation. 
The boundary (matching) conditions at r « a for the 'outer' wave function are related to 
the measurable parameters of the resonance. In the discussed process the BB pairs are 
produced in the P wave. Also beyond the region of strong interaction, i.e. at r > a, there 
is no strong interaction mixing between the "neutral", B°B°, and the "charged", B + B~, 
channels. Thus the 'outer' wave function at r > a is described by the spherical wave with 
L = 1, whose radial part can be written as R(r) = x(r)/r, with a separate function x(r) for 
each of the channels: x„(r) and Xc{ r ), each satisfying at energy E = p 2 /m the corresponding 
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one-dimensional Schrodinger equation 

Xn + ~ 4) Xn = , Xc + + m ~ ~ \) Xc = , (3) 

where the prime denotes derivative over r, and m = tub ~ 5280 MeV" is the mass of either of 
the B mesons (a possible small mass difference between the charged and neutral B mesons 
is completely ignored throughout the present discussion). 

The coupling between the "neutral" and the "charged" channels takes place in the region 
of strong interaction at short distances. At those distances the light quark parts of the mesons 
strongly overlap and become a part of (presumably) quite complicated dynamics of light 
quarks and gluons. Thus in this region it would be inappropriate to continue description in 
terms of individual B mesons. The boundary condition for the 'outer' dynamics at distances 
r > a is however dictated by the isotopic invariance of the strong interaction. Namely, one 
can assume with a rather high degree of accuracy that when the B mesons emerge from the 
region of strong dynamics as individual particles their wave function is an isotopic singlet. 
In other words, the isospin condition for the functions \ n and \ c is that they evolve from 
one and the same function at a certain short distance r = a, i.e. that 

Xc(a) = Xn(a) and x'M) = Xn( a ) , ( 4 ) 

which boundary conditions can be viewed as our formal definition of the short distance 
parameter a. Although there can be a small 'intrinsic' isospin violation also in the region 
of strong interaction, its effect in SR c ^ n , as discussed in Ref.jEl, is noticeably smaller than 
that of the Coulomb interaction, and can be studied as a further adjustment, using the 
approximation of exact isospin symmetry at short distances in eq.(@J) as a starting point. 

As is known from the standard Breit-Wigner description of a resonance scattering [0] 
at energy E near the position E of the resonance, the relevant 'outer' solution of the 
Schrodinger equations (JHJ) for stationary wave functions has the form 

Xn (r) = (A- l7 )6n/nW + (A + H)&;/;(r), 

X c(r) = (A-* 7 )& c / c (r) + (A + z 7 )&:/ c *(r) , (5) 

where A = E — Eq, 'j = T/2 and the complex coefficients 6 n ( c ) are generally functions of the 
energy, which however have no zeros at A = % 7. Finally, each of the functions /„ and f c is 
the solution of the corresponding equation in (jSJ), which contains only the outgoing wave, 
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i.e. at r — > oo they contain only the factor exp(ipr) (while their complex conjugates /*, c ) 
contain only the incoming wave factor exp(— ipr)). 

The function f n {f) specified by this condition is well known for the free motion with 
L = l, 

Ur) = { 1 + ^) ^ ' (6) 

and with this condition for its phase, the phase of the coefficient b n coincides with the 
non- resonant scattering phase 5i at L — 1, 

ex P (2z5 1 ) = £f . (7) 

The corresponding function f c {r) for the motion in the Coulomb potential is also well known 
(see e.g. in Ref. 0), however for our present purpose it would be more convenient to make use 
of the perturbation theory in the Coulomb interaction, rather than to do an expansion of the 
explicit expression. In specifying the phase convention for the function f c a minor technical 
point arises due to the well known fact that its phase at r — > oo contains a slowly varying 
logarithmic Coulomb phase: / c (r) ~ expi[pr + (a/2v) hi2pr + const}. This however can be 
readily resolved by assuming that the Coulomb interaction is cut off at a large distance r = R. 
Then at still larger r the Coulomb phase does not change and can be considered as constant. 
Clearly the physical results, including the discussed here effect in b~R c l n do not depend on 
this infrared cutoff. Implying such regularization, the function / c (r) in eq.([5|) can be chosen 
to exactly coincide (both in phase and in normalization) with / n (r) at asymptotically large 
distances: / c ( r ) = exp(zpr). Thus any difference in phase and magnitude that arises from 
the Coulomb interaction in the "charged" channel is encoded in the coefficient b c . 

In order to find from the wave function of a stationary state (eq.©) the relative rate of 
production of the pairs of charged and neutral B mesons in e + e~ annihilation, i.e. by a source 
localized well inside the region of strong interaction, it is necessary to note that the rate in 
each channel is proportional to the inverse of the norm squared of the coefficient in front of 
the incoming wave. This can be understood by considering the reverse process: annihilation 
of a meson pair into e + e~, in complete analogy with an explanation of the "|?/>(0)| 2 rule" 
for production of bound states. In the reverse process the incoming wave has a fixed flux, 
i.e. a fixed norm of the coefficient in front of exp(— ipr) at large r. Matching this wave to 
the incoming part of the wave function in eq.(jHJ) implies that the corresponding function 
x(r) has to be divided by the coefficient of its incoming wave part. Under the normalization 
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conventions adopted here for the functions /(r) this results in the annihilation rate being 
proportional to the factor |(A + b*\~ 2 , with b equal to b n or b c , depending on the chosen 
initial state for the incoming wave. Clearly, this factor in the rate contains both the standard 
Breit-Wigner resonance curve and the normalization factor \b\~ 2 . Thus the discussed here 
ratio of the yields in the two channels is given by 1 

Rc/n = M ■ ( g ) 

\bc\ 2 

The relation between the coefficients b n and b c is found from the matching conditions 
(jlj). After substituting the wave functions from eq.© into the conditions (jH), the ratio of 
the coefficients (in fact the inverse of that entering eq.(jH}) is found as 

I, _ r: (fn + g e-^W-^j _ g fj, n + r: e-MBw-m^ 

bn fc fc fcfc 

'fc {fn + f: e- 2 ^-- 2 *) - / c * (/: + /;* e - M ™- M >)] , (9) 



2p 

where all the functions and their derivatives are taken at the matching point r = a, 5bw 
is the standard Breit-Wigner resonance phase: exp(2i5sw) = (A — ij)/(A + ij), and the 
non- resonant phase 5± is defined by eq.(J7J). Finally, in the last transition a use is made of the 
fact that the denominator in the intermediate expression is the Wronskian, which is constant 
in r and can thus be found from the asymptotic form of the function f c at large r. 

The equation Q contains no approximation with regards to the Coulomb interaction, 
and can be used down to arbitrarily small values of v, i.e. for arbitrary values of the Coulomb 
parameter a/v, provided that the exact Coulomb function is used for / c (r). However for the 
practical purpose of discussing the Coulomb effects at the T(45) resonance it is sufficient to 
consider only the effect of first order in a. In this order one can write / c (r) = f n { r ) + 4>{ r ) 
with being formally a small perturbation of order a of the wave function. Using this 
expression in eq.® and also assuming an expansion of the ratio b c /b n to the first order in 
a, one readily finds 

6Rc/n = M _ 1 = I l m {eW^+W [0(a) /» - 0'(a) f n (a)}} . (10) 



1 In a somewhat more widely familiar non-resonant situation when point-like particles are produced by a 
point source, the functions x«( r ) an d Xc{r) are the regular at r = solutions of the equations Jjyl. In this 
case one has (for a P wave): \b n /b c \ 2 — \ip' c (0) / ip' n (0)\ 2 , where ip c {r) and ip n (r) are the wave functions for 
corresponding stationary states, having the same relative normalization at infinity. 
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The combination of the functions in this expression can be found directly from the equations 
(J3J) by also a rather standard method. Indeed, the function f n (r) satisfies the first of those 
equations, while f c (r) satisfies the second equation. Multiplying the first equation by f c {r) 
and the second equation by f n { r ), subtracting the results, and expanding in the difference 
between f c and f n , one arrives at the relation 

£ W) f'n(r) - /„(r)] = m - f n {r) . (11) 

Integrating this relation over r from a to infinity, using the fact that under our conventions 
4>(r) — > at r — > oo, and also using the explicit expression (JHJ) for f„,(r), one finds in terms 
of eq.(JTUJ) the formula 



5R c/r ' 



-— Im 

v 



,j2iS bw +2i<5i 



3 2ipr 



i \ dr 
1 + — — 



a: 



A 2 



7 



A 2 + 7 2 



(A cos 25i + B sin 25 x 



r 

2 7 A 
A 2 + 7 2 



(12) 



(5 cos2<5i - A sin2<5i; 



where in the latter expression the coefficients A and B are given (with the oposite sign) by 
respectively the imaginary and the real part of the integral with complex exponent: 



.4 
B 



pa 



1 



1 ) sin 2u 



pa 



2 sin 2m 



u 



u- 



2 cos 2m 
u 

cos 2u 



du 



u 



du 
u 



At small values of the product pa the coefficients A and B have the expansion: 



.4 



7T 

2 



2pa 



3„3\ 



+ 0(p A a 



B 



1 



\n2pa — + 1 + 0(p a 



(13) 



(14) 



3 ■ - v - ' ' " 2 p 2 a 2 

where 7^ = 0.577 ... is the Euler constant. 

One can see from eq. (|12j) that in the limit, where the effects of the strong interaction in 
the wave function of the scattering state are ignored, corresponding to 7 — >• and 5± — > 0, 
the simple estimate (J2J of 5R c ^ n is recovered, assuming production of point-like particles by 
a point source (i.e. also a — > 0). However this limit of a vanishing resonance width is totally 
inadequate for resolving the behavior of 8R c ^ n inside the resonance curve, i.e. at the energy 
scale of 7. In particular, the second term in the final expression in eq. ()12)) . proportional to 
7A/(A 2 + 7 2 ) vanishes at energies far away from the resonance, i.e. at |A| 3> 7, as well as 
at the center of the resonance, i.e. at A = 0. On the other hand, the first term contains 
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the factor (A 2 — 7 2 )/(A 2 + 7 2 ), which changes from +1 away from the resonance to —1 at 
A = 0. Thus bR c l n necessarily has to change sign within the width of the resonance. 

Given that the coefficient B is singular at small momenta, it is instructive to analyze 
the behavior of bR c l n described by the equations (jT2l - IT%j) in the limit of small p (but still 
considering a/v as small for the applicability of the perturbative treatment of the Coulomb 
interaction). At this point one has to recall that the width parameter 7 for a P wave 
resonance has to vanish at small p as p 3 , and that the non-resonant scattering phase S\ also 
vanishes as p 3 . Thus the effect of the p~ 2 singularity in the coefficient B results in a constant 
at small velocity term in the correction: const ■ a, which is small in comparison with the 
dominant part of the correction behaving as a/v. 

It should be emphasized, prior to discussing the behavior at the realistic T(4S) resonance, 
that the expressions (fT2j) - (fTl|) are formally applicable only in the limit of small pa. Indeed, 
in this limit the details of the transition between an isotopically symmetric strong dynamics 
at short disatnces and the Coulomb behavior at the relevant distances of order p~ x are not 
essential, and the parameter a enters only the leading singularity of the coefficient B resulting 
in a sub leading at low energy term in SR c ^ n , which can be studied phenomenologically 2 . 
When the parameter pa cannot be considered as small, the details of the actual behavior of 
the wave functions at the transition distances would generally depend on the shape of the 
transition, and a more elaborate matching at short distances, than the conditions Q at a 
fixed distance, might be required. 

The position of the actual T(4S) resonance corresponds to the momentum po pa 330MeV. 

Thus, most likely, the relevant values of the parameter pa are of order one, i.e. the parameter 

is neither small nor large in the region of interest. Under these circumstances one can 

approach this region of phenomenological interest from the side of small values of pa aiming 

at an at least qualitative description of what kind of behavior should be expected for 5R c ^ n at 

the resonance. In addition, the width of the T(4S), T pa lAMeV, is not small as compared 

to its excitation energy. Combined with the not too small value of po this can generally 

lead to that effects of higher terms of expansion in p 2 of the quantities 7(2?), Si(p), and 

A(p) may become essential. (An attempt at taking into account higher than p 3 terms in 

the width parameter was done in connection with the experimental measurement of V in 

2 It can be noticed that a non-removable dependence on a short-distance parameter is a general property 
of scattering in the states with L ^ 0, and in fact it also appears, although as a logarithmic dependence, in 
the S wave scattering in the presence of a Coulomb interaction|U]. 
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Ref. [10] ■ However the effect of those higher terms turned out to be quite small within the 
experimental accuracy.) 
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Figure 1: The dependence of the ratio R c l n on the excitation energy E = y/s — 2m b in the 
region of the T(4S) resonance (the center position is assumed to be at Eq = 20 MeV) for some 
values of a and 5x(E ): a' 1 = 200 MeV, <J x (Eb) = (solid), a" 1 = 400 MeV, Sx(E ) = 
(dashed), a" 1 = 300 MeV, <$ x (Eb) = 30° (dashdot), and a" 1 = 300 MeV, 6i(E ) = -30° 
(dotted). 

With all the stated reservations about uncertainties involved in a quantitative description 
of the behavior of 5R c ^ n in the region of the T(45) resonance, a qualitative illustration of 
the expected variation of R c / n in the resonance region is provided by the plots in Fig.l. 
The curves in the plots are calculated with various values of a and 5\ under the following 
assumptions: only the leading terms in the expansion of A and B at small pa, explicitly 
shown in eq.(|14|). are retained, the width parameter is parameterized as 7 = (r/2) {p/po) 3 , 
and the non-resonant phase as 8\ = 5i(Eq) {p/po) 3 - (It should be mentioned that the range 
of Si(Eo) from —30° to +30° most likely is unrealistically broad, and is used here for an 
illustration of the effect of the phase under extreme assumptions.) One can clearly see 
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from the shown curves that being undoubtedly different in details, they exhibit quite similar 
qualitative behavior 3 . Thus, as expected on general grounds from eq.(|12|). the very fact of 
a substantial and rapid variation of 5R c ^ n within the resonance width stays robust under 
assumptions about the presently unknown parameters. 

The strong variation of 5R c ^ n at the scale of few MeV near the center of the T(4S") 
resonance can be important in a comparison of the data obtained in different experiments, 
especially at different electron-positron colliders, given the differences in the beam energy 
spread and possible slight differences in the values (and stability) of the central energy 
of the beams, at which the data are collected. The thus far available results ^ 12 E] are 
consistent with each other mainly due to the substantial statistical and systematic errors. 
This agreement however can change if the measurements of the discussed relative yield 
are pursued further with better accuracy, possibly permitting a quantitative study of the 
variation of this yield. As is argued here, the effect of the variation is not small, and the 
detailed behavior is sensitive to properties of the heavy and light quark and hadron dynamics, 
which would be difficult, if possible at all, to study by other means. 

I gratefully acknowledge useful discussions with Ron Poling, Jon Urheim, and Arkady 
Vainshtein of the experimental and theoretical issues related to the problem considered in 
this paper. This work is supported in part by the DOE grant DE-FG02-94ER40823. 
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